Jamming transition in a homogeneous one-dimensional system: the Bus Route Model 
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We present a driven diffusive model which we call the Bus Route Model. The model is defined on a 
one-dimensional lattice, with each lattice site having two binary variables, one of which is conserved 
("buses") and one of which is non-conserved ("passengers"). The buses are driven in a preferred 
direction and are slowed down by the presence of passengers who arrive with rate A. We study 
the model by simulation, heuristic argument and a mean-field theory. All these approaches provide 
strong evidence of a transition between an inhomogeneous "jammed" phase (where the buses bunch 
together) and a homogeneous phase as the bus density is increased. However, we argue that a strict 
phase transition is present only in the limit A — > 0. For small A, we argue that the transition is 
replaced by an abrupt crossover which is exponentially sharp in 1/A. We also study the coarsening 
of gaps between buses in the jammed regime. An alternative interpretation of the model is given 
in which the spaces between "buses" and the buses themselves are interchanged. This describes a 
system of particles whose mobility decreases the longer they have been stationary and could provide 
a model for, say, the flow of a gelling or sticky material along a pipe. 
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I. INTRODUCTION AND MODEL 

Driven diffusive systems |^] have recently attracted 
much attention in the field of non-equilibrium statisti- 
cal mechanics from a fundamental viewpoint, as well as 
in the context of traffic modelling 0, interface growth 
1^ and other applications Q]. 

One particularly interesting feature is the possibility of 
phase ordering and phase transitions in one dimensional 
(Id) systems. To appreciate the significance, one should 
recall that in Id equilibrium models, ordering only occurs 
either in the limit of zero temperature [e.g. kinetic Ising 
models or deterministic Ginsburg-Landau equation) or in 
mean-field-like models. However, in non-equilibrium sys- 
tems it has been demonstrated that ordering may occur 
in models with fully stochastic, local dynamics |^-0. 

The non-equilibrium transitions found so far appear 
to be of three main types (although very recently, novel 
phase separation phenomena have been demonstrated in 
some Id systems P,Hl). First there are boundary induced 
transitions [pO|,p|. These occur on open systems with a 
dynamics that conserves some quantity in the bulk, but 
that allows injection and extraction of the quantity at 
the boundaries. A second class of transition, describing 
roughening of a Id interface, is connected to directed 
percolation and corresponds to a driven diffusive sys- 
tem with non-conserved order parameter |ll| , |l2| . Finally, 
there are transitions induced by defect sites P3|,[l4| or 
particles |15|-p^ or the presence of disorder [19|-|21|j. In 
this class of systems, the presence of the defect causes 
a macroscopic region of high density to form. Analogies 
with Bose Einstein condensation pOl and phase coexis- 
tence of a gas [|l^ have been made. An even simpler 
way to view the phenomenon is as a jamming transition; 
the defect causes a traffic jam to form behind it. In this 



context, however, "jamming" may be a somewhat mis- 
leading term, since in the models just described the in- 
homogeneous "jammed" phase arises at low density. The 
transition is in fact between this phase and a higher den- 
sity "congested" phase which is uniform, but which has 
a mean particle velocity lower than that of the jammed 
phase. Indeed, it appears that a minimum velocity prin- 
ciple applies |jl9|, so that the stable phase is always the 
slowest available at a given mean density. 

In this work we address the question of whether simi- 
lar "jamming" transitions can occur in Id homogeneous 
systems, i.e. systems with periodic boundary conditions 
and without disorder or defects. We introduce a model 
that exhibits a jamming transition in a certain limit (to 
be specified below). For the moment, it is useful to de- 
scribe the model in terms of a commonly experienced 
and universally irritating situation. Gonsider buses mov- 
ing between bus-stops along a bus route. Glearly, the 
ideal situation is that the buses are evenly distributed 
along the route so that each bus picks up roughly the 
same number of passengers. However, owing to some 
fluctuation, it may happen that a bus is delayed and the 
gap to the bus in front of it becomes large. Then, the 
time elapsed since the bus-stops in front of the delayed 
bus have been visited by the previous bus is larger than 
usual and consequently more passengers will be waiting 
at these bus-stops. Therefore the bus becomes delayed 
even further. At the same time, the buses behind catch 
up with the delayed bus and pick up only very few pas- 
sengers since the delayed bus takes them all. Hence a 
"jam" of buses forms. Inspired by this scenario we shall 
formulate a model below, to be referred to as the Bus 
Route Model (BRM). 

We defer the mathematical definition of the BRM un- 
til after we have discussed the general context. Already, 
from the simple picture discussed in the previous para- 



graph, we can identify a conserved variable (the buses) 
and a non-conserved variable (the passengers). The pas- 
sengers are non-conserved since they arrive at the bus- 
stops from outside the system (bus route) . The conserved 
buses arc driven in a preferred direction. Wc may usefully 
think of the non-conserved variable coupling to the con- 
served variable and mediating the jamming transition. 

Associated with any ordering dynamics is the phe- 
nomenon of coarsening p3| where the typical domain 
length of the ordered phase grows indefinitely with time. 
Indeed, coarsening has been studied in ballistic aggrega- 
tion models p3] and disordered driven diffusive systems 
where jamming occurs pO,p3]. Again, a contrast can be 
made with Id equilibrium models where only zero tem- 
perature models or mean-field-like models coarsen. In 
the present model, it is the gaps between the jams that 
coarsen as the jams aggregate; we study this phenomenon 
in the present work. 

In a finite system, the coarsening eventually results 
in one large jam with a single gap in front of it. Re- 
calling that the model system is homogeneous and that 
no bus is preferred over any other, we see that we have a 
spontaneous symmetry breaking where the symmetry be- 
tween buses is broken through one bus being selected to 
head the jam. Symmetry breaking transitions have been 
previously been found in Id open systems with a non- 
conserved variable at the boundaries MM and in a class 
of growth models ||ll|]. The present model provides an 
example in a homogeneous system with a conserved vari- 
able. Related symmetry breaking has also been noted 
in some models in a class of "Backgammon" or "balls 
in boxes" models which are effectively simple generalisa- 
tions of Bose systems 25-27 . However, in these models 



the dynamics is inherently equilibrium and mean-field- 
like whereas the BRM has local dynamics that does not 
satisfy detailed balance. We shall elucidate the connec- 
tion between the two classes of models by showing that a 
mean-field approximation to the BRM results in a model 
that may be solved analytically. The steady state of this 
soluble mean-field model falls into the class of generalised 
Bose systems. 

We now formally define the BRM. The model is defined 
on a Id lattice with periodic boundary conditions. Each 
lattice site is labelled by a number i running from I to 
L. Site i has two binary variables r^ and (j)i associated 
with it. These variables can be described in the following 
terms: 

• If site i is occupied by a bus then r^ = 1; otherwise 
T. =0. 

• If site i has passengers on it then 0j = 1 ; otherwise 

Each site can be thought of as a bus-stop on a bus route. 
A site cannot have both Ti = (j)i = 1 {i.e. it cannot have 
a bus and passengers) . 

There are M buses and L sites in the system and the 
bus density 
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FIG. 1. Schematic illustration of the BRM. 



p = M/L 



(1) 



is a conserved quantity. However, the total number of 
sites with passengers is not conserved. The update rules 
for the system are as follows: 



1. Pick a site i at random. 

2. If Ti = and (pi — then 
A. 



1 with probability 



3. If Ti — 1 and r^+i — 0, define a hopping rate /x as 
follows: 



• ^ = a if (f>^+i = 

• H = P if (t)i+i = 1 



and update r^ 
probability fi. 



n+i 



1 and 
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with 



Thus, a is the hopping rate of a bus onto a site with 
no passengers and /3 is the hopping rate onto a site with 
passengers. The probability that a passenger arrives at 
an empty site is A. When a bus hops onto a site with pas- 
sengers, it removes the passengers. While we have taken 
the passenger variable (f>i to be binary, this does not for- 
bid the presence of more than one passenger at a site; we 
merely require that the extra passengers have no further 
effect on the dynamics. We generally take (3 < a, reflect- 
ing the fact that buses are slowed down by having to pick 
up passengers. We may set a equal to 1 without loss of 
generality and from now on we consider only this case. 
We note that the dynamics is local and does not satisfy 
detailed balance. The model is illustrated schematically 
in Figure 111. 

Although the language of "buses" and "passengers" 
provides an appealing mental picture, the model is in- 
tended to be simple rather than realistic. For example, 
the need for buses to stop to allow passengers to disem- 
bark is ignored. Note however, that an ability for buses 
to overtake each other would have almost no effect. This 
is because, in a jammed situation of the type discussed 
below, the interchange of a fast-moving bus with a slower- 
moving one in front also interchanges their velocities. 

The main purpose of this paper is to provide evidence 
that the model defined above undergoes a jamming tran- 
sition. However, we argue that a strict transition (i.e. 
a singularity in some measured quantity) only occurs as 
A ^ and in the thermodynamic limit. Here we define 
the thermodynamic limit as 



M, L 



with p held fixed. 



(2) 



Our evidence is both numerical and analytical. In Sec- 
tion |l^ we present Monte Carlo simulations which pro- 
vide evidence for the above picture. In the following sec- 
tions we provide various a naly ses to support the picture. 
In particular, in Section II A we present a simple two 
particle approximation which describes the stability of 
jams. In Section III we define a mean-field approxima- 
tion. Within this approximation we can solve analytically 
for the steady state and find it is similar to a generalised 
ideal Bose gas or Backgammon model. This steady state 
can be analysed and exhibits a phase transition (taking 
the form of a condensation transition) only as A —> in 
the thermodynamic limit. We show that the mean-field 
steady state agrees quantitatively with Monte Carlo sim- 
ulations of the BRM, suggesting that for the BRM also 
there is no strict transition for nonzero A. We also study 
numerically the behaviour as A —> in the mean-field 
approximation. 

In Section IV we study the approach to the steady state 
in the jammed phase where we observe coarsening of bus 
clusters. We argue that on an infinite system the size of 
the large gaps between clusters should eventually grow as 
ti/2. We study finite systems numerically. In Section ^ 
we discuss an interpretation of the BRM wherein we con- 
sider the vacant sites between buses to be the moving en- 
tities in the system. Viewed this way, the non-conserved 
variable now describes an internal degree of freedom of 
the moving entities themselves, rather than of the sites 
they visit. We discuss possible physical interpretations of 
this dynamics, including a model of clogging. In Section 



VI, we consider the relevance of the BRM to real bus 



routes and in Section VII, we conclude with a discussion 



of the main points of our work. 



II. SIMULATION RESULTS AND HEURISTIC 
ARGUMENTS 

The model defined above captures an important fea- 
ture of the bus route problem described in the intro- 
duction, namely that once a gap between buses becomes 
large through some fluctuation, the tendency is for the 
gap to become still larger: since buses move more quickly 
in areas with few passengers, a bus which is following 
closely behind another will tend to move faster than one 
which is a long way behind the bus directly in front. This 
is simply because the closer a bus is to the one directly 
ahead of it, the less time passengers will have had to 
arrive. If this tendency for large gaps to grow were to 
prevail, the result on a finite system would be a single 
jam of buses and one large gap. We first argue that in 
the limit of A — > but XL — > cx) this scenario can only 
hold at low enough density of buses, and that a phase 
transition to a homogeneous phase will result as the den- 
sity is raised. 

First consider a system comprising a single large jam. 
In order for this to be a stable object, the velocity (de- 



fined as the average rate of hopping forward) of the lead- 
ing bus must be equal to the velocity of any bus inside 
the jam. Now if XL — > oo, the probability that the site 
immediately in front of the leading bus has a passenger 
on it tends to one. (This is because the rate of passenger 
arrival, multiplied by the time delay between the final 
bus of the jam and the leading bus of the jam crossing 
the same site, is of order XL). Therefore the velocity of 
the leading bus will approach /3. On the other hand, if 
A ^ the probability that a site within the jam has pas- 
sengers tends to zero since these gaps are of finite length. 
Therefore, since the velocity of a bus in the jam will only 
be limited by the presence of neighbouring buses, this 
velocity will be 1 — pjam where pjam is the density of 
buses within the jam[j. Equating the two velocities yields 



Pja 



l — (3. However, the jam must clearly have density 



greater than p, the overall density of the system. Thus 
for jamming to occur we require 



p < Pc where pc = Pjam = 1-/3. 



(3) 



For densities above this critical value, the system will 
be in a congested phase where gaps between buses are 
uniform. 

An equivalent way of obtaining (0) is to compare the 
velocity in the jammed phase, /?, and the velocity in the 
homogeneous phase, 1 — p. The phase with the lowest 
velocity is chosen by the system. This procedure is anal- 
ogous to the thermodynamic procedure of choosing the 
phase with the lowest chemical potential. Though un- 
proven for the present problem, in the case of disordered 
exclusion processes the analogy has been shown to be 
exact I p^ . 

We now present simulations which, for small A, quali- 
tatively support the above picture of a phase transition. 
Figure || shows a space-time plot of the buses in system 
at p = 0.55 for small A (A = 0.02). The buses are dis- 
tributed fairly homogeneously throughout space. There 
are no very large gaps present in the system. Figure y 
shows a space-time plot for p = 0.2 which is less than pc 
of (H) . In this case, starting from a random configuration 
of buses, large gaps quickly open up and small clusters or 
"jams" of buses are readily seen. Gradually, these small 
jams coarsen until, finally, the system comprises a single 
large jam. There is one large inter-bus gap in front of 
the jam whereas the jam itself contains many small gaps. 
The behaviour can be thought of as phase separation into 
regions of nonzero and zero density. 

In order to investigate the effect of varying A, we plot in 
Figure § the steady state average velocity w as a function 
of the density. A system size L = 10000 was chosen 
since using a bigger system did not appreciably affect the 



*This follows from the fact that, when A — > with XL — > oo, 
the situation is equivalent to a model of hopping particles 
with a single slow "defect" particle [[l9[. 




FIG. 2. Space-time plot of bus positions for A = 0.02, 
p = 0.55, (3 = 0.5 and L = 500. There are 10 time-steps 
between each snapshot on the time axis. Initially the buses 
are positioned randomly and there are no passengers. 
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FIG. 3. Space-time plot of bus positions for the same pa- 
rameters as in Figure H except that here, p = 0.2. 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 




A & A fi fi — A — fi — A — ^ 



0.0 0.2 0.4 0.6 

P 



0.8 



1.0 



FIG. 4. Velocity as a function of density for various values 
of A and /3. The simulations were performed with L — 10000. 
The lines are shown to guide the eye. 



results. Let us first consider the results for f3 = 0.5. For 
A — 0.1, the velocity increases smoothly as the density is 
decreased and appears to approach f3 for small density. 
There is no sign of a phase transition. For A = 0.05 a 
similar picture holds although now the curve v{p) is more 
concave. 

A strikingly different picture is obtained for A = 0.02. 
There is an apparent discontinuity in the derivative of 
v{p) at some value of the density p* ~ 0.45. Below p* , 
the data is consistent with v = (3 whereas above p* , v 
decreases almost linearly with increasing density. This 
behaviour is consistent with the simple picture of a jam- 
ming transition discussed above. It is also very similar 
to the velocity-density relationship in an exactly solvable 
model with a single slow particle p9| , p0[ where a jamming 
transition occurs. The graph therefore suggests that a 
phase transition may occur for small A. We argue, how- 
ever, that the phase transition occurs only in the limit 
A — > and that this limit has to be taken in an appro- 
priate way. In order to quantify this we now analyse a 
simple two-particle approximation to the full system. 



A. Two-Particle Approximation 

Consider a system containing M buses. Let us assume 
that there is a jam in the system (i.e. a gap with size 
0{L)). If there is a jamming transition then such a gap 
should become stable in the thermodynamic limit; the 
bus at the head of the jam should not be able to escape 
into the large gap. Let us consider the two buses A and 
B at the head of such a jam as shown in Figure ^. The 
gap in front of A has size kL where k is independent of 
L. The gap in front of B has size x. We now assume 



u(x) 



u(kL) 
(a) Gap of size kL 



where 



FIG. 5. The two buses at the head of a jam in a system 
of size L. 



that we can write the hop rate of either bus as a function 
only of the gap size in front of that bus. To do this 
we write u{x) = f{x) + (3{1 — f{x)) where f{x) is the 
probabihty that there are no passengers on the first site 
of a gap of size x. To estimate f{x) we assume that bus 
A is a random walker hopping with rate w, where v is the 
average velocity of the system. Then the average time 
since A left the site in front of B is x/v. Since the rate of 
arrival of passenger is A we then have f(x) — exp(— Acc/w) 
and 

u{x) — j3 + {\ — P) exp(— Ax/w) for a; > 

w(0) = 0. (4) 

The use of expression (jj) for the hopping rates is in the 
spirit of a mean-field approximation. We have found by 
simulation that the approximation is a good one for j3 
larger than about 0.2, but for small (3 it breaks down. 
In particular, for small A (when j3 is small) we find that 
u{x) decays to /3 much more rapidly than (^ predicts. 
We believe that the reason for this is the failure of (Q) 
to take into account time correlations in the hopping of 
buses - when a bus is updated and fails to hop into an 
unoccupied site, the next time it is updated it should 
hop with probability (3 and not u{x) (because for a bus 
to fail to hop, the site ahead must contain passengers^) . 
Clearly, the effects of neglecting these time correlations 
will be largest for small /?. We shall return to (Q) in 
Section [II where we carry out a conventional mean field 
theory for the many particle problem. For our present 
purposes, we replace w by /? in (Eh since we are interested 
in a jammed situation. 

Using the mean field hopping rate u{x) we may write 
a Langevin equation for the dynamics of the gap size 



x = u{kL) — u{x) + 7]{t) 



(5) 



where 77 (i) is a noise term (say white noise of unit vari- 
ancd'^'h. This can be written in the form 



- — ^{x)+r^{t) 



(6) 



^This is no longer true for 07^!, but time correlations will 
still be present. 

•'•The variance of the noise should strictly depend on j3 but, 
since we are primarily interested in the effect of (small) A on 
the dynamics of the gap, we ignore this dependence. 



$(x) = -(l-/3) 



-XkL/f} 



f3 



-Xx/fj 



(7) 



The gap size x has the dynamics of a particle undergoing 
diffusion in a potential $(a;) for a; > 0. There is a reflect- 
ing boundary at a; = and the potential has a maximum 
at x* — kL. Therefore for x < x* the particle is trapped 
in a well {0 < x < x*) and bound to the reflecting bound- 
ary, i.e. A is bound to the head of the jam. If a; > a;*, 
the particle has escaped from the well and A is no longer 
at the head of the jam. We denote the time to escape 
from the well by r which, to a first approximation P8[ , 
is given by 

T~exp($(2;*)-$(0)). (8) 

In the thermodynamic limit L ^ 00, we have 

"/3(l-/3) 



T ~ exp 



A 



(9) 



which is finite for A > 0, showing that particle A will leave 
the head of the jam in a finite time. This implies that 
a jam is not a stable object and will eventually break 
up. However, when A — *■ 0, particle A becomes bound 
to particle B since r diverges. Therefore, in this limit, 
the jam is a stable object. This simple analysis suggests 
that there is a phase transition (between a jammed and 
a homogeneous phase as the density is raised) only in 
the limit A ^ 0. When A is small but nonzero, r is 
exponentially large in 1/A and it can appear that a jam 
is stable when in fact it has a finite lifetime. 

It is clear that the limits A — > and L ^ cx) do not 
commute since as A —> on a finite system one recov- 
ers a model where all hopping rates are the same. Such 
a model of hopping particles with hard core exclusion 
(known as an asymmetric exclusion process) has been 
well studied and exhibits no phase transition with peri- 
odic boundary conditions pi. In practice, one could take 
the limit L —>■ 00 and then A ^ by choosing 



A - L-"^ with 7 < 1 



(10) 



and taking the thermodynamic linrit L ~> 00 whereupon 
the escape time diverges as 



T '^ exp(aL''') 



(11) 



and the jam is a stable object. 

In order to test this picture against simulations, we plot 
in Figure Ha space-time plot for /? = 0.1 and A = 0.02 for 
a low density. At first glance, it appears that jamming 
just like that seen for /? = 0.5 is taking place. However, 
closer inspection of the individual jams shows that they 
can develop large gaps and, in some cases, divide into 
smaller jams. The steady state for this system is there- 
fore not characterised by a single large jam (as appeared 
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FIG. 6. Space-time plot of bus positions for A = 0.02, 
p = 0.2, /3 = 0.1 and L = 500. There are 20 time-steps 
between each snapshot on the time axis. Initially the buses 
are positioned randomly and there are no passengers. 




FIG. 7. Space-time 
rameters as in Figure 
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lot of bus positions for the same pa- 
except that here, A = 0.002. 



to be the case from simulations for /3 = 0.5), but com- 
prises a number of jams of varying sizes. 

Figure shows a space-time plot for the same value of (3 
(0.1) but for A = 0.002, a factor of 10 smaller. It appears 
that the jamming has been restored ~ the individual jams 
are not seen to break up. It is also interesting to note 
that splitting of jams is not observed for /3 = 0.5 (Figure 
|3|) so the splitting in Figure Ig is a direct consequence of 
having a smaller value of /3. Figure also shows v{p) for 
A = 0.02 and (i = 0.1. Unlike the case /3 = 0.5, v{p) is a 
smooth function with apparently no discontinuity in its 
derivative. This is consistent with the observation that 
the jammed state is unstable on the time scale of the 
simulation (see Figure ^) . 

We therefore conclude from our two-particle argument 
that although the simulation data is indicative of the 
presence of a phase transition, it cannot be taken as firm 
evidence that a strict transition occurs for any nonzero 
A. However, for small A something closely approaching a 
transition is seen. 



III. MEAN-FIELD MODEL 

In this section we discuss a mean-field theory which 
we believe describes the behaviour of the BRM very well. 
Motivated by the mean-field probability (0) of hopping 
into a gap of a given size, derived in Section II A , wc ap- 
proximate the BRM by a model which is exactly solvable 
in the steady state. We term this new model the mean- 
field model (MFM). We analyse this MFM and show that 
it compares quantitatively well with the BRM and, im- 
portantly, that it does not have a phase transition for 
A > 0. However, we find that there is a transition in the 
MFM in the limit A — + so long as the thermodynamic 
limit is taken first. This supports our view that the same 
applies in the BRM itself. 

We define the MFM to be a hopping particle model 
(generalised asymmetric exclusion process ||l|) with the 
hopping probability of the i"^ particle being u(xi), a func- 
tion of Xi, the size of the gap in front of that particle. As 
in the BRM, we consider M particles and L sites in total, 
and periodic boundary conditions apply. The hopping 
probabilities u{x) are given by the mean- field expression 
(^). Note that u{x) depends on v which is the steady 
state mean velocity of the system, (u). It is useful to 
think of V in (Q) as an adjustable parameter in the MFM 
to be determined self-consistently; i.e. v is chosen (for 
fixed p, (3 and A) so that {il)— v. 

While u{x) is given by (Q) for the MFM, it is conve- 
nient to first consider a more general u{x)^ with the only 



constraints being m(0) = (exclusion) and u{x) -^ /3+ as 
X ^ oo. When one considers the dynamics of the gaps 
between particles, it is evident that such a model is an 
example of what is known in the mathematical literature 
as a zero-range process [^,0 which is exactly soluble 
(in the steady-state). We present the solution in detail 



in Appendix |^ and quote the result here. Each config- 
uration can be uniquely identified by a set of gap sizes 
{xi} = {xi, X2, . . . , xm}- The steady-state probability of 
a configuration {xi} is 



Pi{^^}) 



Z{L,M) 



q{xi)q{x2)---q{xM) (12) 



where 



1 

9(2;) = n ~r\ 



(13) 



and 



Z{L, M) 



X! <l{xi)---q{xM)5x,+...+XM,L- 

Xl ,X2,---,XM 



M 



(14) 



Z(L, M) can be viewed as the partition function of a 
generalised non-interacting system of Bose particles. We 
follow Bialas, Burda & Johnston |22l in the consequent 
analysis, starting from the expression for Z{L,M). We 
first discuss the behaviour for general u(x) in the thermo- 
dynamic limit and state the conditions which u{x) must 
satisfy in order for a phase transition to occur. Turning 
to the specific case of the MFM, we show that there is 
a phase transition in the limit A -^ 0. We describe a 
method for analysing finite systems and we compare the 
MFM with simulation results for the BRM. Finally, we 
discuss the behaviour for small, but nonzero, A. 



where 



li^) 



F'{z) 
F{z)- 



(18) 



We may identify z as the fugacity. Each value of the 
fugacity gives a particular value of the density. The par- 
tition function becomes 

Z{p, M) - exp {M log F{z) - [L - M) log(z)} (19) 



with z for a particular p determined by solving (17). 

To obtain quantitative results from this solution, it is 
useful to obtain an expression for p{x)^ the steady state 
probability that a given particle has a gap of size x in 
front of it. For any L and M , p{x) is given by 



p{x) 



q{x) 



Z{L,M) 



q{x) 



Z{L 



y^ q{x2)- ■■q{xM)5x-, 

2,...,XM 

x-l,M -1) 



-\-...-\-xm,L — AI—x 



Z{L,M) 



(20) 



This can be determined for a finite gap size x in the 
thermodynamic limit using the saddle-point expression 
for Z{L,M) given in (|l|). We obtain 



p{x) 



q{x)z^ 

-T{z) 



(21) 



This expression is in terms of z but p may be found using 

(0)- 

The mean particle velocity v in the steady state is 



A. Mean-Field Model in the Thermodynamic Limit 



(") 



E 



u{x)p{x) 



(22) 



Using the integral representation of the Kronecker 
delta 



ds 



2tt\ s"+i 
we write Z{L,M), defined in (p^), as 



Z{p,M) 



ds 
27ris 



Fis) 



-iM 



(15) 



where p = M/L and we have defined the generating func- 
tion 



Fis)=Y,q{x)s^ 



(16) 



x=0 



The integral (^_5h is calculated in the thermodynamic 
limit (defined in (^) using the saddle point method. It 
can be shown that the saddle point, where s = z, is given 
through the expression 



- - 1 = zg{z) 
P 



(17) 



Substituting ( pT| ) into (22) gives the result that v — z. 
This is a relationship which has been found before in this 
kind of system ||T^ . We show below that z is constrained 
to be no greater than (3 and hence, there is an upper 
limit on the velocity. This constraint on z combined with 
(pi| ) implies that in the thermodynamic limit, p{x) is a 
monotonically decreasing function of x and behaves for 
large x as p{x) oc {z/ (3Y . 

We now examine the criteria that u{x) must satisfy 
for a phase transition to occur. This entails analysis of 
the generating function F{z) defined by (16). F{z) con- 
verges for z < /9 and diverges for z > [3 (since we have 



u{x) -^ (3^ 



cx)). Now, g{z = 0) = which cor- 



responds io p — 1. Observe that g{z) is a monotonically 
increasing function for Q < z < (3 which means that p 
decreases monotonically from 1 to 1/(1 -|-/3g(/3)) for z in 
this range. It is clear then that g{l3), if finite, will give a 
critical density, pc, via (p^). Indeed, Bialas, Burda and 
Johnston [g6| show that in this case there is a transition 
from a high density "congested" phase to an inhomoge- 
neous "jammed" phase for p < Pc where a single gap will 
occupy a finite fraction of the sites in the system. The 
critical density is given by 



= lim 



1 



fi- l + zg{z)' 



(23) 



For a transition to occur at nonzero density, we must 
have g{z) (and hence both F{z) and F'{z)) convergent 
as z — > I3~ . 

For convenience, we write u{x) = /3(1 + C{x)). By 
considering the asymptotic behaviour of the product q{x) 
in the large x hmit, one finds that a phase transition 
occurs (i.e. pc is nonzero) if and only if ^(a;) decays 
to zero more slowly than 2/x as a; -^ oo. For the MFM, 
where u{x) is given by (y), we have C,{x) ~ exp(— Aa;) and 
we immediately see that there is no transition for A > 0. 
We now show that, in contrast, in the limit A ^ 0+, a 
transition does occur. 

For the MFM, F{z) is given by 



^w=s i ni 



1 



ay) 



(24) 



where C,{x) = (1//3 — 1) exp(— Aa^/u). To show that a 
transition occurs in the limit A -^ 0^, we must show 
that 



B. Finite Systems 

It is possible to analyse the MFM for finite systems 
since, using (EQ), we have 



L-M 



L-M 



E p(^) = 1 = E 



x=0 



x=t) 



Z{L~x-l,M -l)q(x) 

Z{L, M) 



which gives the following recurrence relation for Z(L, M): 



L-AI 



Z{L, M) = Y^ Z{L-x-l,M ~ l)q{x) 



(27) 



with the initial condition Z(L, 1) = q{L — 1). 

In principle, one can calculate Z{L, M) for any L and 
M . However, numerical precision restricts the practical 
maximum values of L and M to ~ 2000. In practice 
for these finite system calculations, it is not feasible to 
retain the dependence of u{x) on v and so we simply 
take V = 13 va the expression (^ for u{x). This is a good 
approximation since we are principally interested in the 
case where v is close to (3. 



lim lim Flz) 



and 



lim lim F'(z) (25) 

z~^p- A-+0+ 



both converge and also that 



lim lim F(z) = lim lim F'(z) = oo. (26) 

z^/3+ A^0+ z^f3+ A^0+ 

Since for A > both F{z) and F'{z) diverge when z > (3, 
the requirement (Pq ) is satisfied trivially. Note that if the 
limit A — > is taken before the thermodynamic limit, val- 
ues of the fugacity greater than /3 are permitted and we 
recover a model where all particles hop with probability 
1. 

We now show that the expressions in (ESh are finite 
and we calculate the critical density. Clearly, F{z) and 
F'{z) both converge for A > when z < (3 (since the 
geometric part {z/ f3Y dominates the series for large x). 
Since 1 + C,{x) = 1/(3 for A = 0, one can easily calculate 
F{z) — Y^z^ and F'{z) ~ J2^^^~^- Using the saddle 
point equation (|l7|), one finds that the density p is 1 — z. 
Therefore, the critical density is non-zero and is given by 
Pc = 1 — p. We anticipated this result in the discussion 
in Section || in the context of the A — > transition in 
the BRM. Note that the analysis given here does not 
make any predictions about the behaviour for densities 



below Pc. However, the results presented in Sections IIIC 



and HID, together with the discussion in the context of 
the BRM in Section |l[ lead to the conclusion that the 
low density phase comprises a jam in which the leading 
particle hops forward with probability /3 (because the 
site in front of it contains passengers but no bus), and 
all following particles hop forward with probability 1, so 
long as the site in front contains no bus (which holds with 
probability /?). 



C. Comparison with BRM 

We now proceed to compare MFM results with BRM 
simulation data. All the results presented for the MFM 
are based on a calculation of the gap size distribution 
p{x). It is not possible to obtain a closed analytic ex- 
pression for p{x) and so one must perform the calculation 
numerically for a given pair of parameters f3 and A. 

Figure |j shows v{p) for the BRM and the MFM; the 
agreement is quite good. For A = 0.1, we have found 
numerically that in the BRM, system size is not an im- 
portant factor for large enough systems and v{p) for 
L = 1000 is essentially the same as v{p) in the thermo- 
dynamic limit. (The major part of the small discrepancy 
between the MFM data for L = 1000 and L = oo at high 
density can be attributed to the fact that the calculation 
for the finite system is performed with v being replaced 
by P in the expression for u(x).) 

Recall from Figure El that in the BRM, where we had 
L = 10000, for A = 0.02 and (3 = 0.5, v{p) exhibited an 
apparent discontinuity in its derivative, a possible signal 
of a phase transition. For the MFM in the thermody- 
namic limit with the same parameters, v{p) is very simi- 
lar (see Figured) but we know that it is actually a smooth 
function and that there is no transition. This is consis- 
tent with our view that there is in fact no phase transition 
in the BRM for A > 0. Figure g also shows that, for this 
value of A (unlike A = 0.1), there is a marked finite size 
effect at L = 1000 comprising a bump in v{p). The bump 
appears for both the BRM numerics and the MFM cal- 
culation. The reason for the presence of the bump is that 
the size of the large gap in front of a jam decreases as 
L decreases, resulting in the "head" of the jam catching 
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FIG. 8. Velocity as a function of density for the BRM (sim- 
ulation) and the MFM (calculation) for A = 0.1 and A = 0.02. 
I3 is 0.5 and L is 1000. MFM data for the thermodynamic limit 
is also shown. 



up with the "tail" (and hence an average velocity greater 
than P) for densities sufRcently close to pc — I ~ f3- 

Figure O compares the BRM and MFM gap size dis- 
tributions for small systems in the jammed regime. The 
distributions are bimodal with the approximately Gaus- 
sian second peak corresponding to the presence of a single 
large gap in the system. The agreement between simu- 
lation and MFM is again reasonably good. The MFM 
distribution in the thermodynamic limit is also shown 
and one can see that it decays monotonically, but very 
slowly, for large gap sizes. We have found that the posi- 
tion of the peak in the tail of p{x) increases linearly with 
system size as one would expect - the peak corresponds 
to the presence of an "extensive" gap (a gap with size 
(xL). 



D. Small A behaviour of MFM 

We now examine in detail the behaviour of the MFM 
when A is small (but nonzero) . Figure |l^ shows the ef- 
fect on the MFM gap size distribution of decreasing A 
towards zero in a finite system. For A = 0.05, p{x) de- 
creases monotonically from a maximum at x = 0. As A is 
decreased, p{x) becomes bimodal, signalling the presence 
of jams. The second peak becomes more pronounced as 
A is made smaller; jamming is enhanced. 

In order to quantify the behaviour of p{x) as A — > 0, 
we define Pext as the area under the peak in the tail of 
p{x). Then P^xt is the probability of finding an extensive 
gap in the system and Pext x M is the average number 
of extensive gaps in the system. Figure ^ shows plots 
of Pext X M against M for various values of A. As A is 



10" 



10' 



10"' 



- L=1000(BRIVI) 

" L = 500(BRIVI) 

MFM (thermodynamic limit) 

MFM (finite size) 



10° 



10 




0.0 



200.0 400.0 600.0 

Gap size x 



FIG. 9. BRM (simulation) and MFM (calculation) gap size 
distributions on a linear-log scale for small systems in the 
jammed regime {/3 = 0.5, A = 0.02, p = 0.2). The MFM 
distribution in the thermodynamic limit is also shown. 
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FIG. 10. MFM gap size distributions on a linear-log scale 
for a small system showing the effect of decreasing A. The 
other parameters are (3 = 0.5, p — 0.2 and L — 1000. 
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FIG. 11. Pext X M against M for the MFM showing the 
effect of decreasing A. The other parameters are 13 — 0.5 and 
p = 0.1. 

decreased from 0.05, one must go to larger systems to 
observe Pext x M falling below 1. This suggests that 
in the limit A ^ 0, a single extensive gap survives in 
the thermodynamic limit, supporting our claim that a 
condensation (or jamming) transition does occur in this 
limit. 



We showed in Section III A that for the MFM in the 



thermodynamic limit, p{x) is proportional to {v / (3Y for 
large x. Therefore, the typical size of the large gaps in 
the system is the decay constant of p{x) , given by 



e = 



ln(^ 



(28) 



Figure |l^ shows a linear-log plot of S, against 1/A for 
p = 0.3. We see ^ ^ exp(a/A) (where a is a constant) 
for small A so that the typical size of the large gaps is 
very large for A less than about 0.02 and has an essential 
singularity as A ^ 0. Since p{x) is sharply peaked aX x — 
(see Figures || and |l^) , we deduce that for low density 
and small A, a MFM system comprises large clusters of 
buses which are typically a distance ^ apart. Clearly, this 
can only apply if L ^ ^; if this is not the case, then the 
gap size distribution must be bimodal as seen in Figure 

Figure |l^ shows the effect on v{p) of decreasing A 
towards zero in the thermodynamic limit. As A is 
decreased, a "corner" becomes apparent at a density 
slightly below 0.5. This behaviour is also observed in 
the BRM (see Figure 0) . While we know that there is no 
singularity in v{p) for A > 0, one can see that as A ^ 0, 
v{p) — > /3 for p < 0.5 and v[p) — > 1 — p for p > 0.5. We 
now examine the sharpness of the crossover from the low 
density jammed regime to the high density homogeneous 
regime for small A. 
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FIG. 12. Linear-log plot of the typical large gap size ^ 
against 1/A for the MFM in the thermodynamic limit. The 
density is fixed at p = 0.3 and j3 is 0.5. The lines are shown 
to guide the eye. 
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FIG. 13. Velocity as a function of density for the MFM in 
the thermodynamic limit showing the effect of decreasing A. 
(3 is 0.5. 
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FIG. 14. Linear-log plot of Kmax against 1/A for the MFM 
in the thermodynamic limit. /3 is 0.5. The lines are shown to 
guide the eye. 



To estimate the sharpness of the crossover for a given 
set of parameters, we calculate K^ax, which we define 
to be the maximum value of |w"(p)|. Figure 04 shows 
a linear-log plot of K^ax against 1/A. For A less than 
about 0.02, we see that Kmax goes as exp(6/A), where b is 
a constant. Therefore, although a strict phase transition 
occurs only in the limit A — > 0, the crossover is exponen- 
tially sharp in 1/A for small A; this may be compared 
with the typical large gap size ^ discussed above which is 
also exponentially large in 1/A. For practical purposes, 
the crossover may be indistinguishable from a phase tran- 
sition as is already the case for A = 0.02 (see Figures 
and |l^). Note that, as mentioned previously, the ther- 
modynamic and A ^ limits do not commute: for a 
finite system the behaviour at A = is trivial (there are 
no passengers in the system). This means that in a large 
but not infinite system, the sharp crossover will, as A is 
reduced, resemble a phase transition most strongly for 
some nonzero A, X*{L), and then fade away for A ^ A*. 
(Heuristically, we expect X*L to be of order unity.) 



IV. COARSENING 

In this section we discuss the approach to the steady 
state in the jamming regime. We have already seen in 
Figure H that, in finite systems, coarsening of bus clus- 
ters (or equivalently the gaps between clusters) occurs for 
small A and low bus density. Here, by relating the model 
to a reaction-diffusion process, we argue that the typi- 
cal size of the large gaps in the system should eventually 
grow as t^' ^ and we provide numerical evidence for this. 
We also study numerically the relaxation of the average 
particle velocity v(t) and find that it decays as a power 



law with an exponent close to —1. 

The gap size distribution in the jamming regime can be 
considered as a superposition of the small gaps inside bus 
clusters and the large gaps between these clusters. We 
define r(t) to be the mean size of the large gaps. Recall 
first from Section [II that in a mean-field approximation 
we expect the probability that a particle hops into a gap 
of size r to be 



u{r) ^ /3+{l- /3) exp(-Ar/w) 



(29) 



If Xr is sufficiently large, the motion of a typical bus 
cluster therefore becomes uncorrelated with the motion 
of the cluster ahead. When two clusters come suffi- 
ciently close together, they coalesce and form a single 
cluster. The larger r{t) becomes, the more this coales- 
cence is diffusion-dominated {i.e. correlations in the mo- 
tion of clusters are reduced) and the more it resembles 
a driven A + A ^ A reaction-diffusion process [B2l with 
bus clusters taking the place of A particles. In the driven 
A + A ^ A proccssF], the characteristic length scale (of 
which r{t) is an example) grows as t^^'^. We therefore 
expect t^/^ growth in the BRM when Xr{t) is sufficiently 
large, i.e. at sufficiently late times. 

We have performed simulations on large systems (L — 
5 X 10^) to investigate the approach to the steady state. 
In these simulations, inspection of the complete gap size 
distribution shows that, after a short time, there is a very 
deep minimum at a gap size of about 20; we observe a 
negligible number of gaps of this size. Therefore, in the 
results presented below we have defined a large gap as 
one with size greater than 20 and r{t) is then the mean 
size of gaps greater than 20. In discussing the coarsening 
in a disordered driven diffusive model, Krug and Ferrari 
pO{ use the variance of the complete gap size distribu- 
tion as their measure of the typical length scale. We have 
also studied this quantity and we find that its time evolu- 
tion is entirely consistent with that of r{t), as one would 
expect in a scaling regime. 

Figure 15 shows log-log plots of r(t) obtained from sim- 
ulation. Close inspection of the r{t) curves reveals that 
they are straight lines at early times indicating power law 
growth with exponent ~ 0.3. By early times we mean the 
interval 3000 < t < 10000 for the two larger values of A 
and 10000 < t < 30000 for A = 0.005. However, at later 
times we see different behaviour as the curves become 
somewhat concave, particularly for the two larger values 
of A. 

Since we anticipate that t^'^ growth will occur only 
after r{t) has become sufficiently large for bus clusters to 
be uncorrelated, we expect r{t) to be of the form 



r(i) = ro(A)+Ati/2 



(30) 



^For the A + A ^ A process, the presence of a preferred di- 
rection does not change the scaUng from that of the undriven 
system ]3l[. 
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FIG. 15. Log-log plots of the average large gap size r{t) 
for small values of A. Simulations had L — 500000, p — 0.2 
and (3 = 0.2. Results for A = 0.01 and 0.005 were averaged 
over 15 independent runs and for A — 0.008, 8 runs. At 
t — the buses were positioned at random and there were 
no passengers. The dashed lines (f^'^ and t"'^) are shown to 
guide the eye. 
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FIG. 16. Log-log plots of r{t) — ro for small values of A. 
The data is the same as that presented in Figure [151 The 
dashed line (f^ ) is shown to guide the eye. 



at late times. The constant A should depend only very 
weakly on A if the coarsening process is truly diffusion- 
dominated. However, rg is a measure of how far apart bus 
clusters must become to be uncorrelated for a particular 
value of A, and so we expect it to vary as ~ 1/A from 

Figure |l^ shows log-log plots of r{t) — ro for the same 
data presented in Figure nS. The parameter tq was 
estimated by fitting the function ( |3(]| ) to the data for 
t > 40000. For the two larger values of A we see that 
the curves collapse at late times as expected. The be- 
haviour at these late times is certainly consistent with 
^1/2 growth. However, the results for A = 0.005 do not 
fit this picture. We believe that this is because the t^'"^ 
growth regime has not yet been reached for this value 
of A. It appears that as A becomes smaller, the regime 
characterised by approximate i"'^ growth increases in du- 
ration. Our results suggest that the BRM may exhibit a 
crossover from power law growth of r{t) with exponent 
~ 0.3 to growth with exponent 1/2. We remark that 
the MFM exhibits similar, but not identical, coarsening 
behaviour. 

We now discuss the behaviour of the average velocity 
v{t) as the steady state is approached. We believe that 
the true steady-state velocity in the jammed regime of the 
BRM is very close to, but slightly smaller than, /3. Figure 
[l7| shows the decay of v{t) towards its steady-state value. 
The data, although noisy, shows clearly that v{t) relaxes 
towards f3 = 0.2 as a power law with an exponent close 
to —1. The data for later times, although very noisy and 
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FIG. 17. Log-log plots of v(t) — P for small values of A. 
Simulations are the same as those for the data shown in Fig- 
ures ha and h6l The dashed line (i~^) is shown to guide the 
eye. The data has been rescaled vertically for clarity. 
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not shown, is consistent with t~^. We have no evidence 
for a crossover in the relaxation of v{t) (in contrast to 
r{t)) but, owing to the quahty of our data, neither can 
we rule out this possibility. 

It is interesting to note that we observe coarsening in a 
system which does not strictly phase separate. We have 
already argued that for small A the thermodynamic sin- 
gularities in the limit A — > are replaced by behaviour 
that is smooth but has crossovers that are exponentially 
sharp in 1/A. Apparently the A ^ phase transition 
likewise shows up in the dynamics where the system ap- 
pears to phase separate, but only up to a finite (but very 
long) time so that a truly inhomogeneous state is never 
reach ed. Our study of the MFM for small A in Section 
[IIP suggests, however, that the final "homogeneous" 
state comprises large clusters of buses separated by gaps 
having a typical size of order exp(l/A); hence t^^'^ growth 
of the typical large gap size will occur up to some time 
exponentially large in 1/A (for a sufficiently large sys- 
tem). Perhaps not surprisingly, this timescale is of the 
same order of magnitude as our estimate of the timescale 
on which individual clusters are stable, derived within 



the two-particle approximation in Section [I A 



V. DUAL MODEL: CLOGGING 

The above completes our study of the BRM. However, 
as mentioned in Section |, our motivation for studying 
the model is mainly connected with its interesting generic 
behaviour (a driven system with one conserved and one 
non-conserved variable) rather than its applicability to 
public transport. In this section we present an alterna- 
tive interpretation of the model which further reveals its 
generic behaviour. 

One can interpret the BRM in a different way by noting 
that each time a bus hops to the right, the site (hole) 
that the bus hops into moves to the left. By considering 
the "holes" (henceforth we call them particles) to be the 
moving entities in the model and the "buses" to be empty 
sites, one defines a dual model. 

A feature of the dual model is that the non-conserved 
variable can now be thought of as being "attached" to the 
moving particles rather than fixed sites as in the origi- 
nal interpretation. The non-conserved variable is /i^ , the 
speed of particle i, which can either be fast (/i^ = 1) or 
slow {fii = j3). If the dynamics of this dual model is to 
be exactly the same as that of the BRM, then a particle 
should attempt to hop to the left when the site to its left 
is updated. However, a more natural dynamics is to only 
choose particles for update. The full update rules in the 
latter case are: 

1. Pick a particle i at random. 

2. If /ii = 1 then ^i ^ (3 with probability A. 

3. If there is no particle on the site to the left of par- 
ticle i, then it hops to the left with probability ^i. 



4. If particle i hops, then /i; — > 1. 

This model describes a system of particles each of which 
can exist in two states of mobility. A particle has prob- 
ability A of switching to the less mobile state each time- 
step for which it remains stationary. When it finally does 
move, it is restored to the more mobile state. 

While the dual model is not identical to the BRM, 
we have checked that the numerical behaviour is indeed 
very similar. The density of particles is 1 — p = p' and 
we note that a gap in the BRM is equivalent to a cluster 
of adjacent particles in the dual model. Thus, the mean- 
field theories for the two models are equivalent, so long 
as the hopping rate into a gap of size x in the MFM is 
interpreted as the hopping rate of a particle leaving the 
left edge of a cluster of size x in the dual mean-field the- 
ory. Jamming becomes a high density phenomenon here, 
characterised by the presence of large clusters of parti- 
cles. This restores to the word "jamming" a meaning 
closer to that used in everyday life. 

We define the average speed of particles as v' . Since 
in the BRM the magnitude of the bus current must be 
equal to the magnitude of the hole current, we have 



P 



1 



P 



v{p) = 



1-p' 



v{l^p'). 



(31) 



The dual model can be thought of as describing stop- 
start traffic flow with the particles representing cars. The 
longer a car is at rest (usually because a car in front is 
blocking it), the more likely it is that the driver will be 
slow to react when it is possible to move again. Now A 
is the parameter determining the strength of the "slow- 
to-start" behaviour of cars. This is similar to several 
"slow-to-start" cellular automaton traffic models studied 
recently fS^]. In the study of traffic fiow, the so-called 
"fundamental diagram" (the current v' p' as a function of 
density p') is commonly examined. Figure |l8| shows the 
fundamental diagram for the dual model (simulation and 
MFM) for different values of A. One can see that for high 
densities and A very small, the current is significantly less 
than that for A = 0. It follows from the behaviour as 
A — > that in the model, an infinitesimal probability for 
cars to become slow-to-start results in a macroscopically 
large decrease in the current at high density. 

A different interpretation of the dual model provides a 
simple picture of another familiar kind of jamming which 
might be called "clogging" . Consider particles suspended 
in a fluid which is pumped along a narrow pipe. Imagine 
that the particles, which are comparable in size to the 
pipe width (so they cannot pass each other), have some 
tendency to stick to the walls of the pipe, but that this 
process takes time to occur and therefore can only hap- 
pen if a particle remains stationary with respect to the 
pipe wall for a significant time. A stuck particle will re- 
main stuck if there is another such particle in front of it. 
If not, the stuck particle will detach from the wall and 
move on, but only after a delay. This offers an interpre- 
tation of the dual model in which the attachment rate of 
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FIG. 18. Fundamental diagram for moving holes model for 
different values of A. For all simulation data, L = 10000 and 
MFM results are in the thermodynamic limit. The uppermost 
solid curve is the MFM result in the thermodynamic limit 
followed by the limit A — > 0. The dashed curve is the exact 
result when A is set equal to zero before the thermodynamic 
limit is taken. The latter two curves are identical for p < 0.5. 



a stationary particle is A and the rate at which a stuck 
particle will detach and move on is /3. (The parameter a 
is set by the flow rate of the fluid.) In the limit A — > 0, 
a phase transition then arises from a homogeneous phase 
in which the particles move quickly, to a jammed phase 
of inhomogeneous, slow-moving particles. This transition 
occurs when the particle density is raised above 1 — pc, 
where pc was defined earlier for the BRM. 

More generally, such a model could be taken as a highly 
simplified discrete description of any fluid that has a ten- 
dency to clog. There are many fluids which will solidify 
when at rest but remain in a fluid state if kept mov- 
ing rapidly enough - examples include colloid/polymer 
mixtures |34], clay gels (used as drilling muds in the oil 
industry) and, in some circumstances, blood. A similar 
description might apply to particulate suspensions (say 
in a horizontal pipe) which, if not kept moving, will settle 
under gravity into a relatively immobile deposit. 



VI. RELEVANCE TO REAL BUSES 

While the BRM is not designed to model a real bus 
route with any great accuracy, it is worth commenting 
on its possible relevance to the bus route problem. To do 
this, we must relate our model parameters to those of a 
real (in our case, Scottish) bus route. L is the number of 
stops on the circular route and the bus density p is the 
number of buses per stop. We expect p to be quite low 
(perhaps 0.1), and hence to be in the regime where jam- 





(a) 



(b) 



FIG. 19. Space-time plots of bus positions for p = 0.2, 
13 = 0.5 and L = 100. In (a), A = 0.1 and in (b), A = 0.02. 
There is 1 time-step between each snapshot on the time axis. 
The system is in the steady state. 



ming could potentially occur. The parameter a, which 
we take to be 1, is inversely proportional to the aver- 
age time a bus takes to travel from one stop to the next 
(with no stopping for passengers). The parameters /? and 
A are defined relative to a. The amount by which buses 
are delayed owing to their having to pick up passengers 
is reflected in (3. Our choice of /? = 0.5 is, we believe, 
a reasonable figure for a city bus route - a bus picking 
up passengers at every stop will progress about half as 
quickly as one which has no passengers to collect. We in- 
terpret A as the probability that one or more passengers 
arrive at a typical bus-stop in the time it takes a bus to 
go from one stop to another (roughly a minute, say). Of 
course, there is no such thing as a typical bus-stop on 
a city bus route but we believe that A should often lie 
between 0.1 and 1, depending on location, time of day 
and other factors. 

In our examination of the BRM up to this point, we 
have primarily studied the limit of large system size L. 
For real bus routes, however, we expect L, the number of 
bus-stops, to be in the region of 50 to 100. Figure |l^ (a) 
shows the positions of 10 buses in a system of 100 bus- 
stops for two hundred timesteps in the steady state, with 
A = 0.1. Transient clusters of buses are clearly visible in 
the system, a scenario familiar to regular users of the bus 
route. A similar space-time plot is shown for A = 0.02 in 
Figure [lS| (b). As was the case for larger systems, we see a 
single large cluster of buses which is relatively stable. Our 
model suggests, therefore, that such catastrophic clusters 
of buses will not form if the arrival rate of passengers is 
sufficiently large, but that at some intermediate arrival 
rates, a given bus is quite likely to be part of a small 
cluster. Of course, as the arrival rate is increased still 
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further, the BRM predicts an increasingly disordered bus 
route with httle clustering. 

Finally, we wish to comment on the implications of 
a "hail-and-ride" bus system, where there are no fixed 
bus-stops and passengers may hail a bus at any point 
on the route. In the BRM, this corresponds to a system 
with a large number of stops (large L) and very small 
values of both A and /3. Hence, one expects the buses to 
cluster very strongly, much as in Figure m (with a much 
lower bus density, however). Compounding this prob- 
lem, with /3 very small the leading bus in a cluster moves 
extremely slowly. Thus, the use of a "hail-and-ride" sys- 
tem, while on the surface more convenient for passengers, 
would have dire consequences for the efficiency of the bus 
service. However, it is worth pointing out that on routes 
used by very few passengers, A might be so small that 
one would observe essentially A = behaviour, i.e. a 
homogeneous system with buses being scarcely delayed. 
This is because the probability that a site has passengers 
on it cannot be greater than XL; if XL ^ 1, no bus is 
significantly delayed by having to pick up passengers. 



VII. DISCUSSION 

In this work we have provided strong evidence that 
the BRM undergoes a jamming transition as a function 
of the density p of buses in the limit where A, the rate 
of arrival of passengers, tends to zero. This provides an 
example of a homogeneous system with local, stochastic 
dynamics that exhibits a jamming transition in one di- 
mension. Although it remains an open problem to find 
an exact solution of the model, we have shown that the 
steady state of a mean-field approximation is solvable. 
This approximation captures the essence of how the or- 
dering into jams occurs - the density of passengers at a 
site provides information about the time at which a bus 
last visited the bus-stop. From the elapsed time it is in- 
ferred how far away the next bus is along the route. Thus 
the passengers mediate an effective long range interaction 
between the buses. The nature of the mean-field approx- 
imation is to replace this "induced" interaction (which is 
subject to stochastic variation) with a deterministic one. 

The transition also has strong analogies with the jam- 
ming transition induced by a single defect particle or dis- 
order [O-pOJ. In that case it has been shown that the 
transition is reminiscent of Bose condensation |Q. In 
the present model an interesting point is that there is 
no defect bus present; all buses are equal. Instead, the 
system spontaneously selects a bus behind which a jam 
forms and the waiting passengers condense in front of 
this bus. As mentioned in Section |, it would make very 
little difference if overtaking of buses were allowed in the 
model, since its effect would be merely to interchange the 
leading two buses of a jam. The new leading bus would 
then proceed as slowly as its predecessor. This makes 
the physics different from the defect mediated case where 



unhindered overtaking would prevent a jam ever arising. 
We have shown, at least in the two body approximation 
of Section [I A , that the time for the lead bus to escape 
from its jam diverges strongly with system size. This 
may be compared with the "flip time" in another model 
exhibiting spontaneous symmetry breaking 0]. 

The bus route model has an interesting dual model, de- 
scribed in Section 0. This considers the spaces between 
buses to be moving entities. If the moving entities are 
interpreted as cars, the dual model is a particular type of 
slow-to-start traffic model [pSJ . If the moving entities are 
particles suspended in a fluid, the dual model is a model 
of "clogging" in the transport of sticky particles, or a 
gelling fluid, down a pipe. In either interpretation, there 
is (in the limit of small A) a phase transition between a 
homogeneous and an inhomogeneous phase; however, the 
inhomogeneous (jammed) phase now arises at high den- 
sity. This contrasts to the BRM itself and related models 
|13|-p0[ of defect-mediated jamming, for which jams arise 
when the vehicle density is too low. 

The dual model is, in our view, interesting because 
it describes the jamming of particles whose mobility de- 
pends on an internal dynamical degree of freedom. This 
is represented by the nonconserved variable ^. In the 
dual model itself, the nonconserved variable keeps an 
internal record of how long it is since the particle last 
moved. More generally however, other models of this 
type could entail a nonconserved internal variable rep- 
resenting, say, the orientation of rodlike particles. (Such 
particles would be much more likely to jam in some orien- 
tations than others.) We have shown that the existence 
of an additional, nonconserved degree of freedom is, at 
least in a specified limit (A -^ 0), enough to cause a 
symmetry-breaking phase transition to a jammed state, 
in a homogeneous one dimensional driven system. 

We now comment on our results for nonzero A. Al- 
though we argue that a transition only occurs as A — > 0, 
we have shown that a strong vestige of the transition re- 
mains for small values of A. In fact, in the mean-field the- 
ory of the BRM, the crossover between the two regimes 
becomes exponentially sharp in 1/A. Therefore, in prac- 
tice the crossover becomes very difficult to distinguish 
from a strict thermodynamic transition when 1/A is of 
order 100 or more. The A — > transition also strongly 
influences the dynamical behaviour of the system when A 
is small but non-zero: we observed coarsening behaviour 
(presumably transient) as would usually be associated 
with systems that do strictly phase separate. All this 
highlights the fact that that care is required to unambigu- 
ously identify phase transitions, as opposed to crossover 
phenomena. Indeed, it is possible that closely related 
crossover phenomena occur in cellular automata models 
of traffic ||,|6|. 

Finally, we show in Figure ^ the average time that a 
passenger is required to wait for a bus. The system size 
is 10000, unphysical in terms of real bus routes but rep- 
resentative of the thermodynamic limit. For A = 0.1, we 
see that, as one would expect, the average waiting time 
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FIG. 20. Plot of the mean passenger waiting time as a func- 
tion of bus density for two different values of A. The simula- 
tion parameters are L — 10000 and /3 = 0.5. 



increases smoothly as the density of buses is reduced. 
However, for A = 0.02, it increases very sharply at an 
intermediate value of the density, and for low densities, 
the bus service becomes highly inefficient. This figure 
serves as yet another cruel reminder of the vagaries of 
the bus route to those whose lives are unfettered by the 
ownership of a motor car. 
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APPENDIX A: STEADY STATE SOLUTION OF 
MEAN FIELD MODEL 

The MFM comprises M particles hopping on a Irf ring 
of L lattice sites. The probability of hopping into a gap 
of size X is u{x) which must obey u(0) ~ (hard-core 
exclusion). The update is random sequential so that the 
probability of a given particle being updated at each step 
is ^. 

A configuration C is completely specified by a sequence 
of gap sizes {xi} where Xi is the size of the gap in front 
of the i*^ particle. 

In the steady state, the master equation for the system 



J2 W{C -^ C')P{C) = Y^ W{C' ^ C)P{C') (Al) 



C'TtC 



C'i^C 



where P(C) is the steady state probability of configura- 
tion C and W(C -^ C) is the probability of an elementary 
transition from C to C 



First, consider the RHS of dAl]). Define C = 
{xi, . . . , Xi, . . . , xm\ and Cj = {xi, . . . , Xj^i - 1, Xj + 
1, . . . , xai}- Then the only configurations from which C 
can be obtained by an elementary transition are the Cj, 
with the constraint that Xj-i > 0. We therefore have 



WiC, ^ C) = jjuix, + 1) 



and the RHS of (|A1|) becomes 

l^J2u{x, + i)p{c,)e{x,.,) 



(A2) 



(A3) 



where 9{x) is the usual Heaviside (step) function. 

Now consider the LHS of (|Ai|). Define Ck = 
{xi, . . . , Xk-i + l,Xk - I,- ■ ■ , xm}- The Ck are the only 
configurations which can be obtained from C by an el- 
ementary transition, with the constraint that Xk > 0. 
Therefore we have 



W{C ^ Ck) = -^uixk) 



and the LHS of ( |AlD becomes 

l-J2u{xk)P{c)e{xk). 



Equating (^) and ( |A^ ) gives 

J2uMPic)e{xk) = ^w(x, + i)p(c,)0(x,_i 



(A4) 



(A5) 



(A6) 



Since periodic bounda ry c onditions apply, the indices on 
the LHS and RHS of (A6) can be matched up to give 



J2u{xj-i)P{C)eixj-i) = ^^i(x, + l)P(C,)0(a:j-i). 



(A7) 

To solve this equation, we assume a product form for 
P. We write 



Pi{^^}) 



1 



Z(L, M) 



q{xi) ...q{xi) ...q{xM) (A8) 



where Z{L, M) is a normalisation. Substituting this 
product form into (A7) gives 



Y.\\{q{xMe{x,^M3=^ 



(A9) 



where 



A, =.(.,_!) -^^i:^i±ll#^V^iM:fi±l). (Aio) 

q{xj-i)q{xj) 
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To solve this for q{xj), make the ansatz that 



for a;,_i > 0. This has the solution 



^(^) = n 4t 



(A12) 



which, together with (|A8[ ), gives -P({a;j}) as required. 
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